Creation of Clusters via a Thermal Quench 



Yossi Farjoun 

Department of Mathematics, MIT 



Summary. The nucleation and growth of clusters in a progressively cooled vapor is 
studied. The chemical-potential of the vapor increases, resulting in a rapidly increas- 
ing nucleation rate. The growth of the newly created clusters depletes monomers, 
and counters the increase in chemical-potential. Eventually, the chemical potential 
reaches a maximum and begins to decrease. Shortly thereafter the nucleation of new 
clusters effectively ceases. Assuming a slow quench rate, asymptotic methods are 
used to convert the non-linear advection equation of the cluster-size distribution 
into a fourth-order differential equation, which is solved numerically. The distribu- 
tion of cluster-sizes that emerges from this creation era of the quench process, and 
the total amount of clusters generated are found. 



1.1 Introduction 

While studying a simplified model of nucleation the temperature is often as- 
sumed to be held constant (see, e.g., the review papers by Wu [5] and Oxtoby 
[S]). However, time dependence of the temperature can play an important 
role in the aggregation process. In this paper we study the effect that a ther- 
mal quench has on aggregation. Together with J. Neu, we previously studied 
[3] a similar problem with constant temperature. For a less terse literature 
summary, the reader is referred to the review articles cited above. 

Our system comprises a dilute, condensable vapor in a carrier gas. Start- 
ing from a temperature corresponding to zero chemical-potential, the system 
is cooled uniformly at a constant rate. We assume a quench rate that is slow 
relative to the molecular timescale. This implies that the maximal chemical- 
potential is small and therefore the nucleation rate and cluster growth can 
be assumed to follow the Zeldovich formula [7] and the Becker-Doring (BD) 
equations [T] respectively. We use the BD equations and not diffusion lim- 
ited growth, as per Lifshitz-Slyozov [4], since the final size of the clusters is 
relatively small. To find the value of the chemical-potential we use conserva- 
tion of particles and the Clasius-Clapyron relation together with the ideal-gas 
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law. To simplify the model and calculation, we assume that the latent heat of 
evaporation is large relative to the thermal energy, fc^T. 

Initially, there are very few clusters and so the chemical-potential increases 
as the temperature drops (through the resulting reduction of the equilibrium 
monomer density). As the chemical-potential increases, so does the nucleation 
rate. Eventually, enough clusters have been created that their growth causes 
a large enough combined drain on the monomer density so that the increase 
in chemical-potential is stopped. After this, the nucleation rate drops quickly 
and approaches zero. 

The paper is organized as follows: In Section[L2]we present a short deriva- 
tion of the aggregation model we use. In Section 11.31 we solve the resulting 
non-linear advection PDE using an asymptotic approximation, the method of 
characteristics, and a numerical solver. 

1.1.1 Assumptions 

Throughout the paper we make the following assumptions. 

1. Nucleation occurs at the Zeldovich rate. 

2. The clusters that form are small, and grow according to BD dynamics. 

3. The process is spatially uniform. 

4. The only interaction between clusters is via the chemical-potential. 

5. The temperature is not affected from the condensation of clusters. 

6. The carrier gas and the condensable gas are ideal gasses. 

1.2 The Model 

We briefly derive the constituent equations of the quenching process. We as- 
sume basic familiarity with standard nucleation, BD growth, and the monomer 
conservation argument that leads to the determination of the chemical- 
potential from the distribution of cluster sizes. For these we follow the no- 
tation in Wu's review article First, we derive the relationship between the 
chemical-potential and the undercooling. 

1.2.1 Chemical-potential and the Undercooling 

The chemical-potential rj is the free energy of a monomer in condensed liquid 
relative to that in vapor: 

r/^fcsTlog-. (1.1) 

Ce 

Here, ks is the Bolzmann constant, T is the temperature, c is the ambient 
monomer concentration, and Ce is the equilibrium monomer concentration. 

The dependence of Ce on the temperature can be found through the 
Clausius-Clyperon relation and the ideal-gas law as in 
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Co T 

where A is the latent-energy of condensation (per monomer.) In (|1.2p . Tq and 
Co are the initial equilibrium temperature and monomer-concentration. That 
is, at concentration co and temperature To the vapor phase is in equilibrium 
with the liquid phase. Thus, setting T = To gives Ce = co. 

We assume that A ^ ksTp, and define e = , . The non-dimensional 
versions of c, T, and 77 are: 

c==— , T = —, V^TT^- (1-3) 
Co Jo Kb Jo 

Equations p.ip and p.2p now are as follows 



77 = Tlogf, 5, = ie---. (1.4) 



c 

Ce 

The undercooling is defined as r = 

Using the above definition of tau we find an asymptotic approximation of 
Ce for small e: 

Ce = e"^ + O(eT), (1.5) 
and the leading order approximation of fj: 

Tywlogc-fT, for e ^ 1. (1.6) 

1.2.2 The Zeldovich nucleation rate and the growth rate of clusters 

Two important components of our model are the rate at which new clusters 
come into existence, and the rate at which existing clusters grow. The BD 
equations of growth specify that the size n of a cluster follows the growth 
"law" for clusters much larger than the critical size: 

h = ujfjn^ . (1-7) 

Where rj is as in (II. ip , and w is a "escape rate" constant so that uj-n?/^ is the 
rate at which monomers leave the cluster. 

The Zeldovich formula give the nucleation rate of new clusters: 

j=ucodeJ-^e (1.8) 

Where a is the "surface energy" constant of a cluster, 

surface energy 82 , , 

Using the definition of r and p.Sp we have 
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1.2.3 Growth dynamics via advection PDE 

We are interested in finding the evolution of the density of cluster sizes. Let 
r{n, t) be the density of clusters of size n at time t (also referred to as the 
"distribution of cluster-sizes"). From equation (|1.7p we derive an advection 
PDE for the distribution of clusters in the space of their size, 

dfT + ujfidn (^n^rj — 0, in n > 0. (l-H) 
The flux of clusters ujfjnir asymptotes to the Zeldovich rate as n — > 0, thus 
Lofjnir ^ j = ujcqa -^e ^ ^ as n ^ 0. (1-12) 

V DTT 



1.2.4 Determination of the chemical-potential 

The chemical-potential can be inferred from the distribution of cluster sizes 
and the initial concentration cq. 

/•oo 

c+ nr{n, t) dn = Co- (1-13) 

"'0 

For small values of r/ we have from (|1.6p that c = cqc''"'^. Thus 



e^' + nr{n, t) dn = 1, where f — — . (1-14) 
Jo Co 



1.3 The Mathematical Problem 



We now drop all tildes and refer only to non-dimensional variables. The math- 
ematical problem is therefore, 

dtr + rjdninir) = 0, (1-15) 




(1.16) 



e''-^ + / nrdr = 1. (1.17) 
Jo 

Here, time t is non-diniensionalized with the scaling unit 1/lu. 

In the current paper we consider an undercooling which increases linearly 
with time t: 

T^f2{t + tQ). (1.18) 

The parameter J7 is externally specified, and to is a time-lag needed so that 
the "interesting" behavior happens near < = 0. 
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1.3.1 Asymptotic Solution of the Creation Era 

The equations needed to find the relevant scales of the creation era are mostly 
straightforward dominant balances of equation (|1.15l - [07)) . There is one that 
is not: The change in chemical-potential must be such that the reduction in 
the nucleation rate is comparable to the nucleation rate itself. This implies 
that the change in chemical-potential is small (relative to rj itself) and we use 
dr] = 77(0) — ?7(t) to follow the change in the chemical potential. To save space, 
we omit the derivation and proceed to the resulting scales. They are 

While to and i7 are connected by il^tQE^ — . We introduced E, a measure 
of nucleation rate: E = ^/^e ^Co")^ Using these scales results in the 

V 67r ^ 

following system for r{n, t) and 77: 

dtr + dn{nir) = Q, (1.19) 
n^r ^ e'^'' as n ^ 0, (1.20) 

6ri{t) -t+ nrdr = 0. (1.21) 
Jo 

Since n^r is constant along the level curves of — i, we can write 

r(7i,t) = n-ie''''(*-3"^). (1.22) 

Substituting this into (|1.2ip yields an integral equation for S-q, and using the 
change of variable t' = t — Sna we arrive at 

Sr){t)~t + J* (^j^^ e^'J^^'Ut' ^ 0. (1.23) 

From here we derive a fourth-order ordinary differential equation (ODE) for 
6r], by differentiating 4 times: 

2 

,577 = --e'^'', where 5?7(t) ~ t, as t ^ -cx). (1.24) 
9 

To solve this ODE (numerically), we start with Sr]{ti) = ti, Srj{ti) = 1, Sri{ti) = 0, 
Si){U) = 0, and t, = -10. 

By integrating e'''' we get the (scaled) density of clusters p that were 
formed, p « 5.1. The resulting density is shown in Fig. 11.11 
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Fig. 1.1. The cluster size distribution r{n,t) for t ranging from to 9 and, inset, 
the nucleation rate as a function of (shifted and scaled) time. 



1.4 Conclusions 

Thermal quench is a standard trigger for nucleation. However, normally the 
quench process itself is ignored and only the outcome (i.e., a super-saturated 
monomer solution) is considered. In this paper we have shown that during 
a slow quench, enough clusters nucleate so that no more nucleate after the 
quench. Presumably, the clusters that have nucleated simply grow after the 
quench and eventually coarsen. 
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